Phase ordering and shape deformation of two-phase membranes 
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Within a coupled-field Ginzburg-Landau model we study 
analytically phase separation and accompanying shape defor- 
mation on a two-phase elcistic membrane in simple geometries 
such as cylinders, spheres and tori. Using an exact periodic 
domain wall solution we solve for the shape and phase or- 
dering field, and estimate the degree of deformation of the 
membrane. The results are pertinent to a preferential phase 
separation in regions of differing curvature on a variety of 
vesicles. 



Amphiphilic molecules assemble in aqueous media to 
form bilayers, which close to form vesicles at low con- 
centrations. Bilayers and vesicles serve as models for 
membranes and cells for studying simple physical proper- 
ties such as shape deformations, elasticity and transport. 
They show an amazing variety of shapes, which have been 
described by treating the membrane as a homogeneous 
elastic sheet with area and volume constraints . How- 
ever, recent experimental observations have recognized 
that internal degrees of freedom can crucially influence 
the shapes. An example is the transition from a nor- 
mal biconcave shape of discocytes to a crenated shape 
of echinocytes of a human red blood cell. Such transfor- 
mations can be induced by an asymmetric adsorption of 
certain drugs, i.e., a local asymmetry in the composition 
plays an important role in this crenated shape 

As molecules are free to move in the plane of the mem- 
brane, lateral phase separation is constantly observed in 
lipid membranes. A single component membrane under 
certain conditions can exhibit regions rich in tilted and 
non-tilted phases, respectively, while a two-component 
membrane can exhibit phase separation of both different 
components and tilted vs. non-tilted phases [Q. Phase 
separation plays a central role in the stabilization of vesi- 
cles and in the fission of small vesicles after budding . 
Although experimental studies that clearly relate phase 
separation to local shape deformation Q| are scarce, a 
number of phenomenological and numerical investiga- 
tions have shown that a coupling of the local curvature to 
the local composition of amphiphilcs can result in shape 
deformation [Q-^ and budding . 

The numerical studies have considered mainly a gen- 
eral fluctuating vesicle and therefore the central role of 
the coupling between the phase separation and accompa- 
nying shape deformation process has been difficult to de- 
cipher. Our work is thus motivated by a desire to study, 
by analytical means where possible, phase separation on 
the simplest of geometries such as cylinders, spheres and 



tori. While microtubules are abundant in biology, de- 
formable spherical and toroidal vesicles have also been 
observed We seek to provide a systematic formal 

description of the equilibrium solution. Our aim is to 
gain insight into the part played by curvature and to ex- 
tract the salient ingredients that affect phase separation 
and shape changes. We estimate the degree of deforma- 
tion from the coupling strength between the composition 
and curvature fields and the elastic rigidity. 

We represent a membrane as a surface embedded in 
three dimensions parameterized hy q = {91,(72}, for its 
thickness is usually several orders of magnitude smaller 
than its size. Our approach is to study phase separation 
on a subset of surfaces (orthogonal curvilinear manifolds 
]l2[) which have either an axis of translation or rotation. 
For this special class of surfaces we have recently derived 
some simplifying analytical results regarding phase sep- 
aration on rigid curved surfaces [|l3| . Here we apply the 
analysis to deformable surfaces. 

The total free energy of the membrane is F = i^i -|- 
P2 + P3, with the bending elastic energy H : 
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where dA = ^Jgd^q is the area element with y/g the deter- 
minant of the metric tensor gij, h the mean curvature, 
and /iQ the spontaneous mean curvature-the preferred 
curvature of the relaxed vesicle. The mean curvature is 
h = /i((7i, (72) = hi + h2, where hi and /12 are the principal 
curvatures. For an arbitrary surface embedded in three 
dimensions, if qi and 92 are orthogonal coordinates, the 
metric tensor has gij = for i ^ j and ^/g = ^giig22- 

We treat phase separation within the usual Ginzburg- 
Landau free energy framework: 
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where (j) is either the relative concentration of the two 
phase components A and B of the membrane: <?!> = (A — 
B)/{A + B), or the concentration of a diffusing external 
chemical, as in the case of echinocytosis of red blood 
cells [§,0. Here ^ is the characteristic length, which 
determines interface width. 

We use the bilinear form of coupling between the phase 
density and local curvature p4[ , an interaction energy 
found on phenomenological grounds: 
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where A is the strength of the couphng. 

More reaUstic considerations should take into account 
area and volume constraints for vesicle membranes. A 
change in area would increase the surface energy and a 
change in volume would increase the osmotic pressure. 
Hence an additional term 



F4, = XA + PV 
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where A is the surface tension and P the osmotic pressure, 
should be included in the free energy F. The second term 
in the expression utilizes the divergence theorem in 3D. 
The surface tension A is a constant and does not enter the 
variational calculations. The relation between the local 
radius \r\ and the mean curvature h can be highly compli- 
cated depending on the geometry, which will render the 
free energy too intractable for our purposes. Therefore, 
we consider only P — hereafter. 

In the case where the order parameter is conserved, as 
we consider no exchange of molecules of the membrane 
with its environment, an appropriate equation for the 
time evolution of the phase ordering field (j) is: 
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where prime denotes differentiation with respect to (p. 
Similarly, for the mean curvature: 

§ = r^vis^ = r^vL Mh - ho) + Acj^] , 

where Fi and T2 are the diffusion coefficients for (j) and h 
fields, respectively, and Vlg = -j^ {9'^ Vd -£t) is the 
Laplace-Beltrami operator. As the details of the double- 
well potential are not important |ll5|| , we use a simple 
(/)^ potential to describe the phase separation dynamics: 

y(<^) = f04_|^2^ («,/?> 0). 

The numerical study of the dynamics of evolution will 
be reported elsewhere Here we focus on the an- 

alytic equilibrium solutions. At equilibrium, the Euler- 
Lagrange (EL) equations for and h fields are SF/S(j) = 
and SF/6h = 0, respectively. The EL equations are 
nonlinear and usually do not have an exact closed-form 
solution. In order to obtain analytical results, we con- 
sider special symmetries to reduce the problem to a 
quasi-one dimensional one. If qi is the axis of symme- 
try, = \d(t)/dq2\. We then define a new variable 
Ti as dri = a/ (722/511 dq2. With this new variable. 

The equilibrium condition, derived from the EL equa- 
tion for h is 
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i.e., at equilibrium the local mean curvature of the mem- 
brane is linearly proportional to the local (j). This linear 
relationship explains why phase separated regions have 



local curvature {A/K)(ka result that appeared in the nu- 
merical study of Ref. ||9[. Thus, we can eliminate h from 
the free energy. It follows that the EL equation for the 
free energy with respect to (j) becomes 



Here Ve is the new effective potential: 
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which depends only on (j). The coefficient of the cf)'^ term 
is renormalized and the effective potential becomes an 
asymmetric double-well due to the linear coupling. 

Twice integrating Eq. (^) we obtain a general periodic 
domain-wall lattice solution: 
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where dr = dri^gn = ^g22dq2 is the arc length vari- 
able; a, 5, c and d are real roots of — Vb — 0, i.e. 
Ve - Vo = .gii(0 - a)((/) - h){(t) - c){(t) ~ d) with a > 
b > (j) > c > d; Vo, To arc two constants of integration. 



^ = ^^yA/a^/2/{a — c)(b — d) is the rescaled characteris- 
tic length scale, k is the modulus of the Jacobian elliptic 
function sn(T, k). The shape of the periodic solution de- 
pends on the modulus k, which in turn depends on all the 
parameters of the model and the boundary condition Vq. 
The value of k ranges between and 1. For fc = 0, the Ja- 
cobian elliptic function reduces to a sinusoidal function. 
For fc = 1, sn(T, fc) changes to a kink solution, tanh(r), 
which is no longer periodic and thus is only allowed on 
an open geometry. For a closed geometry, e.g. a circle, a 
periodic solution is required and the number of periods 
for a fixed perimeter L depends on the value of fc [|l7| : the 
periods allowed should satisfy L/C^ — AmK where K{k) 
is the complete elliptic integral of the first kind and m 
is an arbitrary integer. With the linear relationship be- 
tween (j) and ft., Eq. (||), we also obtain the expression for 
the mean curvature /i(t) = ho — -^(^(r). 

Transforming t to space coordinate 52 and exploit- 
ing the fact that qi is the axis of symmetry, we obtain 
4>{liil2) and h{qi,q2). From h{qi,q2), we can then use 
the relationship between h and r, to obtain r[qi, (72)7 the 
shape in real space coordinates. 

For illustration purposes, we now carry out this formu- 
lation on a cylinder with rotational symmetry r = r{z). 
The metric tensor has geg = r^, gzz = 1+r' and 1702 = 0. 

With 9 as the axis of symmetry, we define the new 
variable ri as 



dri 



dz\ 
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The free energy for such a cyUnder is then 



(9) 



in which Ve{4>) is the same as Eq. (^. Applying the 
EL equation for h, we obtain the same hnear relation 
between h and ^ as in Eq. (||) in terms of Ti. Replacing 
h in the free energy, F becomes a function of (j) only and 
its EL equation with respect to (j) is 
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This equation yields the same solution for 
as a function of the arc variable r: 
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and hence the solution of hir) based on the linear rela- 
tionship between and h. 

In order to convert the and h results to phase dis- 
tribution and deformation on a cylinder, i.e., to obtain 
(j){z, 6) and r{z, 0), we need to change the variable r back 
to the original variable z as follows. We replace r — r{z) 
by p(t), such that 



The mean curvature h{z) in terms of p{t) is 

h{T) = I / p + prr / \/~^~~p^ ■ 



(11) 



A numerical integration of this equation provides p(t) 
and thus z ~ J dry^l — pr^ , which in turn gives r{z) — 
p(r(z)), the deformation along the z axis of a cylinder. 

Figure la shows a typical plot for the equilibrium 0, 
h and r as a function of z, obtained with the following 
parameters: a — 4, f3 — 2, ho = 5, £, = 0.2, k — 0.02, 
and A = 0.03. Figure lb is the corresponding axially 
deformed cyhnder, whose shades of gray correspond to 
the amplitude of the order-parameter field. Applying the 
same formulation to a cylinder with an axial symmetry 
r — r{9), we obtain a cylinder with deformation occurring 
only along the cross section. Figure 2a shows a series of 
deformed circles with periods 3 to 6. A radially deformed 
cylinder, as shown in Fig. 2b, is formed by translating 
the deformed circle along the z axis. Figure 2b uses the 
same parameters as in Fig. 1. 

The degree of deformation, defined as the ratio of 
the maximum to the minimum radii, D = rmax/fmim 
is a quantity that can be measured by reflection in- 
terference contrast microscopy [ p^ or atomic force mi- 
croscopy We estimate the degree of deformation 
of an axially symmetric cylinder, for which r = l/H 
and thus D = {Hq — ^c)/{Ho — ^d), as a function of 
K and A. The value of k can be either obtained by using 
micropipette methods pO|, or estimated from molecular 



dynamics simulations ||2^. Although difficult to mea- 
sure directly in experiments, A may be derived from first 
principle molecular dynamics simulations as in pl[ | by 
using Eq. (||). Figure 3a shows that the deformation D 
has a quadratic form in 1/k. Figure 3b indicates again 
a quadratic form dependence of the deformation on the 
coupling constant A. This can be understood by expand- 
ing D in the small deformation limit, when A <^ k. In 
this limit, c — d and cd remain approximately indepen- 
dent of A and k, thus D oc Ci + C2A/K + CsA^/k^ with 
ci, C2 and C3 constants. Indeed, if we vary both A and k 
such that A/k remains constant, the deformed circles are 
virtually identical (figure not shown). 

We also applied this framework to other deformable 
geometries and symmetries, e.g. a sphere and a torus 
with appropriate g^ . Figures 4a, b and 4c, d show results 
for spheres and tori with axial and azimuthal symmetry, 
respectively. The parameters used are the same as in 
Fig. 1. We observe preferential phase separation similar 
to the case of cylinders. 

To summarize, we have developed a theoretical frame- 
work in which we obtained exact analytical solutions 
for the equilibrium phase distributions and membrane 
shapes, for a series of geometries including cylinders, 
spheres and tori. This framework allows for an esti- 
mate of the degree of deformation from the coupling 
strength and the elastic rigidity of the membrane. Since 
fluid properties are essential for modeling cells and mem- 
branes, our model augmented with a coupling to hydro- 
dynamics will enable the study of realistic biological cells. 

The order-parameter field considered above need not 
be a scalar density or relative concentration field. If we 
choose magnetization (M) or polarization (P) instead of 
(j), we would then expect periodic stripes of magnetic and 
polarization domains in regions of different curvature. 
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FIG. 1. Phase separation and deformation on a radiaUy 
symmetric circular cylinder, (a) Plots of order parameter (p, 
curvature h and radius r as a function of z. (b) Equilibrium 
shape of the deformed cylinder; shades of gray correspond to 
the order-parameter (j). 
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FIG. 2. Phase separation and deformation on an axially 
symmetric circular cylinder, (a) Cross sections of deformed 
cylinder with 3, 4, 5 and 6 modes, respectively, (b) Equilib- 
rium shape of a deformed cylinder with mode 6; shades of 
gray correspond to the order parameter (j). 





FIG. 3. Deformation of the axially symmetric cylinder as 
a function of: (a) Elastic rigidity k for fixed A — 0.02; circles 
are numerical data and the solid line is fit Inset shows the 
deformed cross sections of the cylinder, the inner-most curve 
corresponding to smallest k value, (b) Coupling constant A 
for fixed k = 0.01; circles are numerical data and the solid 
line is a fit to a quadratic form. Inset shows the deformed 
cross sections with fixed perimeter, the inner-most curve cor- 
responds to the largest A value. 
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FIG. 4. Phase separation and deformation on spheres and 
tori, (a) Axially symmetric sphere, (b) Azimuthally sym- 
metric sphere, (c) Axially symmetric torus, (d) Azimuthally 
symmetric torus. Shades of gray correspond to difi'erent mag- 
nitudes of the order parameter 0. 
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